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0^ ■ The tunneling splitting in biaxial ferrimagnetic particles at excited states 

o. 

with an explicit calculation of the prefactor of exponent is obtained in terms 
of periodic instantons which are responsible for tunneling at excited states 

s 

i ' and is shown as a function of magnetic field applied along an arbitrary direc- 

tion in the plane of hard and medium axes. Using complex time path- integral 

O ■ 

we demonstrate the oscillation of tunnel splitting with respect to the mag- 
^ . nitude and the direction of the magnetic field due to the quantum phase 

interference of two tunneling paths of opposite windings . The oscillation is 
gradually smeared and in the end the tunnel splitting monotonously increases 
with the magnitude of the magnetic field when the direction of the magnetic 
field tends to the medium axis. The oscillation behavior is similar to the 
recent experimental observation with Fes molecular clusters. A candidate of 
possible experiments to observe the effect of quantum phase interference in 
the ferrimagnetic particles is proposed. 
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I. INTRODUCTION 

The macroscopic quantum phenomenon in spin system at low temperature has attracted 
considerable attention both theoretically and experimentally for more than a decade |3]-§J. 
The magnetization vector in a single domain ferromagnetic(FM) grain and the Neel vec- 
tor in a single domain antiferromagnetic(AFM) grain can tunnel from a metastable state 
to a stable one, which is called the macroscopic quantum tunneling(MQT), or display a 
coherent oscillation between two degenerate states, which results in the superposition of 
macroscopically distinguishable (classically degenerate) states ( the understanding of which 
is a long-standing problem in quantum mechanics) and is called macroscopic quantum coher- 
ence(MQC). The geometrical phase (known as the Berry phase) interference plays a crucial 
role in the MQC. The quenching of MQC can be interpreted by the quantum interference 
between tunneling paths of opposite windings which possess a phase with obvious geometric 
meaning 0-0]. The quenching of MQC for half-integer spin has been shown physically to 
be related to Kramers' degeneracy, however, the effect of geometric phase interference is far 
richer than that. For example, when the external magnetic field is applied along the hard 
anisotropy axis, a new quenching of MQC occurs and is not related to Kramers' degeneracy 
since the external magnetic field breaks the time reversal symmetry || . The Zeeman energy 
of the biaxial spin particle associated with the external magnetic field produces an additional 
geometric phase of tunnel paths which leads to the quantum interference, and the tunnel 
splitting therefore oscillates with respect to the magnetic field. The oscillations of the level 
splitting for the ferromagnetic particles have been verified by the experiment with molec- 
ular clusters Fes which at low temperature behave like a ferromagnetic particle ||. The 
experimental observation of the oscillation of tunnel splitting has triggered off more detailed 
investigations along this direction [|iO|-|T2"|j . Since the tunneling rate in AFM particles is much 



higher than that in FM particles of the same volume |H| the AFM particles are expected to 
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be a better candidate for the observation of macroscopic quantum phenomena than the FM 
particles. The quantum tunneling of the Neel vector in AFM particles has been well studied 
in terms of the idealized sublattice-model in which the external magnetic field does not 
play a role since the net magnetic moment vanishes. The biaxial AFM particles with a small 
noncompensation of sublattices or in other words biaxial ferrimagnetic particles have to be 
considered in order to obtain the effect of the external magnetic field on the tunnel splitting. 
The oscillation of tunnel splitting at ground state of the biaxial ferrimagnetic particles was 
predicted recently with the magnetic field applied along the hard axis [15|. In the present 
paper we investigate the effect of quantum phase interference at excited states for a biaxial 
ferrimagnetic particle in the external magnetic field applied along an arbitrary direction in 
the plane of hard and medium axis. Since the effect of geometric phase interference has 
been observed in the experiment of Fe§ molecular clusters with the magnetic field along an 
arbitrary direction, the present generalization to the ferrimagnetic particles is not only of 
theoretical but also of practical interests. At ground state one only considers the paths of 
imaginary time under barrier. The extension to excited states is highly nontrivial. Paths of 
complex time have to be taken into account since a path at excited states also approaches 
the region of potential well and therefore is of real time. 



II. EFFECTIVE LAGRANGIAN OF A BIAXIAL FERRIMAGNETIC PARTICLE 

IN A MAGNETIC FIELD 

We consider a biaxial AFM particle of two collinear FM sublattices with a small non- 
compensation. Assuming that the particle possesses a X easy axis and XOY easy plane , 
and the magnetic field h is applied along an arbitrary direction in the plane of the hard axis 
(Z axis) and medium axis(Y axis), the Hamiltonian operator of the AFM particle has the 
form 

H=J2 (kjf + k ^f - lKS z a - ih y Sl) + JS X ■ S 2 , (1) 

0=1,2 
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where k±,k u > are the anisotropy constants, J is the exchange constant, 7 is the gyro- 
magnetic ratio, and the spin operators in two sublattices Si and S 2 obey the usual com- 



mutation relation 



°0! °b 



iheijkSabSj; k = x,y, z; a,b = 1, 2). In order to obtain the 
Lagrangian of the system, we begin with the matrix element of the evolution operator in 
spin coherent-state representation by means of the spin coherent state path integrals 



(N f \e-* 6T ' H \Ni) 



'M-l 



n 



fc=i 



M 



n(-v,| ( - !t " /fi iA', ,) 



k=\ 



(2) 



Here we define \N) = |ni)|n 2 ), \N M ) = \N f ) = \nij)\n 2 j), \N ) = \Ni) = \n lti ) \n 2ji ), 
tf — ti = 2T and e = 2T/M. The spin coherent state is defined as 



n a )=e ie «°«\S a ,S a ),(a = l,2) 



(3) 



where n a = (sin# a cos0 a , sin# a sin</> a , cos# a ) is the unit vector, O a =sm(f> a S^—cos4> a S^ and 
\S at S a ) is the reference spin eigenstate. The measure is defined by 



25 + 1 

dn(N k )= ] [ — " 8m.6 a ,kd0 a ,ku< >„.!,■ ■ 



a=l,2 



4tt 



(4) 



Evaluating the path integral on the right hand side of the Eq.(2) we obtain in the large S 
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{N f \e-* 6T '*\N i )= / J] D[9 a ]D[<j> a ]ex V 



0=1,2 



h 



(Lo + Ljdt 



(5) 



with 



L = S a <p a (cos 9 a — 1) — JSiS 2 [sin 9% sin 9 2 cos (0i — <j> 2 ) + cos 61 cos 9 2 ] 



(6) 



0=1,2 



Li = - ^ (fci^a cos 2 a + fc,,^ sin 2 a sin 2 a - jh z S a cos a - 7/i y 5" a sin 9 a sin O ) , (7) 



0=1,2 



where L + Li denotes the Lagrangian. Sin ce spins Si and S 2 in two sublattices are almost 
antiparallel, we may replace 9 2 and (f> 2 by 9 2 = n — 9\ — e e and (fi 2 = n + 4>i + where e e 
and e denote small fluctuations. Working out the fluctuation integrations over e g and e 
the transition amplitude Eq.(5) reduces to 



{Nf\e- 2i6T ' h \N i )= D[9]D[^]exp 



h 



Ldt 



(8) 
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M 1 + M 2 1 M ■ XL (A2 ^12 ■ 2n 

< > - — < > cos 9 + — - ( : < > sin H 
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V 



where V i 



QK ± (cos0-Mh z /2K A 



flK u sin 2 (sin - Mh y /2K U sin i 



(9) 



and 



(9x, 0i) has been replaced by (0, 0). M a = ~fhS a /Q (a = 1, 2), M = •yh (S x - S 2 ) /fi with fi 
being the volume of the AFM particle and x± = 7 2 /^- = 2k ± S 2 /Q and K„ = 2k tl S 2 /Q 
(setting Si = S 2 = S except in the term containing Si — S 2 ) denote the transverse and the 
longitudinal anisotropy constants, respectively. 

We assume a very strong transverse anisotropy, i.e., K± 3> K u . For this case, the Neel 
vector is forced to lie on a cone of angle 26$. Where cos#o = Mh z /2K± = 5h z /h c (5 = 
K u /K±, h c = 2K n /M ). Introducing the fluctuation variable rj such that = 0q+ r] 
and considering Kj_ ^> K tl we have V(0,<f>) = QK± sin 2 o r] 2 + QK U sin 2 # o (sin0 — b) 2 {b = 
Mh y /2K n sin9o = /zsina /5^/(h c /5) 2 - h\ (sin 2 9 = 1 - (5h z /h c ) 2 ), where a is the angle 
between the magnetic field and Z axis) and thus the Eq.(9) is written as 



— — - + — 0H cos 6»o - A„sm o (sm0-&) 

2il_|_7 z 7 Z / 7 7 



?f - if j_ sin 2 # ( 77 + 



M0 



2KxT1 sin 6*o 



(10) 



Carrying out the integral over rj we obtain 



(iV.le- 2 ^/^,) = / 7J[0]exp ( -\ P L eff dr 



(11) 



where 



(12) 



is the effective Euclidean Lagrangian. t = it and j3 = iT . I = I a + If where If = 
QM 2 / (2 , y 2 K±) and I a = fix± sm2 ^o/7 2 are the effective FM and AFM moments of inertia 
17fl, respectively. V (0) = OiT,, sin 2 ^(sin — fe) 2 is the effective potential and 9 = h(So — d) 



(So = S\ + S2 and d = h z /h = hcosa/h with ho = h/jlf). The second term in the 
Eq.(12) , i.e., has no effect on the classical equation of motion, however, it leads to a 
path dependent phase in Euclidean action. When h y =0, V(4>) = K u Qsm 2 9 sin 2 (ft possesses 
the form of the sin-Gordon potential and the directions with 9 = 9 , <fi = and tt are 
two equilibrium orientations of the Neel vector (Fig. 1(b)) around which the small oscillation 
frequency of the Neel vector is seen to be co> = ^2K U Q sin 2 6q/I. The quantum tunneling 
of the Neel vector through two paths of opposite windings results in the quantum phase 
interference. When h y 7^ 0, the potential V((j>) = QK U sin 2 #o(sin — b) 2 has an asymmetric 
twin-barrier (Fig. 2 (a)), and the net magnetic moment of the uncompensated sublattices in 
the applied magnetic field shifts the equilibrium orientations of the Neel vector to <fi = <p + 
and tt — + (0 + = arcsin 6)(Fig.2(b)) around which the small oscillation frequency of the 
Neel vector is modified as uo = cugVl ~ b 2 . The quantum tunneling of the Neel vector through 
two different barriers leads to the quantum phase interference. 



III. QUANTUM PHASE INTERFERENCE AS H y =0 

When the external magnetic field is applied along the hard axis(Z axis), the effective 
potential is V((p) = QK U sin 2 9q sin 2 <fi . The quantum tunneling at finite energy E is dom- 



inated by the periodic instantons |18|]. From the Euclidean Lagrangian (12), the equation 
of motion of the pseudoparticles moving in the classically forbidden region in the barrier is 
seen to be 

£(g'-VM~* (13) 

The Neel vector may rotate by tunneling through potential barriers from one orientation(0 = 
0) to another (0 = tt) along clockwise path and anticlockwise path (Fig.l). The instantons 
satisfying periodic boundary condition are found to be 

4>f = ±— ± arcsin [kisn(ujQT)] (14) 
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where "-" denotes the clockwise path and "+" denotes the anticlockwise path(see Fig.l), 
sn(o;or) is the Jacobian elliptic function with modulus 



E 



(15) 



The two trajectories of instantons <j)f are shown in Fig. 1(a). The Euclidean actions evaluated 
along the trajectories of periodic instantons are 



Sf = W e + 2Ep + i6l 



W, 



2 V dr 



- vm 



dr 



LO 



Of = [ Q^dr = ±0(tt - 2 arcsin k[) 
J-p dr 



(16) 



(17) 



(18) 



where K,(ki), E(k\) are the complete elliptic integrals of the first and the second kinds, 
respectively, k' 2 = 1 — k\ = E/£lK n sin 2 9$. To investigate the quantum tunneling and 
related quantum phase interference at excited states, we begin with the instanton induced 
transition amplitude 



m,n 



(19) 



Pe is the operator of projection onto the subspace of fixed energy |jT9]| . \E E ) and \E E ) are 
two excited states lying on different sides of the barrier. From Eq.(19) the tunnel splitting 
is written as 



AE 







d(f) f d(j)iip E ((f)f)ijE(<j)i)G((f)f } (3; <f> u -(3) 



G = J V[<j>) exp ^~ J P L eff dr 



(20) 



(21) 



When the quantum phase interference of tunneling through clockwise and anticlockwise 
paths is taken into account, the Eq.(20) is written as 



exn (^) 
AE ^L±\l+ + I-\, 



(22) 



± 



^d<t>tr E {^)M<t>t)G{<t>% (3; <f>f, -(3) = exp v , 



I = J d<j>fd<f>iil;%(<i>f)il)E(<i>i)G(<i>f,P; 0;, 



o- 



(23) 



(24) 



G= / £>[0]exp(-- / L e// rfr 



(25) 



(26) 



Jo is independent of tunnel directions. The phase independent tunneling kernel G is now 
evaluated with the help of the periodic instantons. Following the procedure of the periodic 
instanton-calculation in Refs.(20) and (21) a general formula for Eq.(24) is found to be 



Jo~2/3exp 



2E/3\ 



h J 



hwo 



exp — 



We 
h 



(27) 



To investigate the quantum phase interference at excited state, we have to consider additional 
phases coming from the real-time paths in the potential well between — > <f>i and 2 — > n 
(0 — > —0i and —02 — ► — 7r)- Thus tunnel splitting Eq.(22) is rewritten as 



where 



«P(f) 
^P(f) 7 



li exp I — — + exp 



ft 



exp 



h 



+ exp 



h 

h 



(28) 



^[0,±^i]U 



[±</» 2> ±7r] 



q=26 arcsin k 1 , 



(29) 



(30) 



W± = ±W- 



E - 2V(<f>) 
2 J[0,±4>i]u[±4>2,±M y/E- V(4>) 



d(f>. 



(31) 



It is obvious that W+ = W~ . Substituting Eqs.(27), (29) and (30) into the Eq.(28), we 
obtain the tunnel splitting 

AE = 7^ ex P \^s(An)\ (32) 



/c(jfei) " V n 

where A = Sq — d. The tunnel splitting AE is a function of the external magnetic field and 
energy. 



For low lying excited states(fc' 1 = \/E/VLK n sin 2 9q « 1) in which we are interested, 
the energy E may be replaced by the harmonic oscillator approximated eigenvalues E m = 
(m+ 7;)uJoh. Expanding the complete elliptic integrals K{k\) and E{k\) as power series of k' 
and taking note of limit JC(k[ — ► 0) — >|, we obtain the tunnel splitting of the mth excited 
state, 

AE m =^^AE \ cos(Att)|, (33) 



nil 



where 



AE =^(8B)ieM-B) (34) 



with B = 4K U Q sin 2 9q/1tujq which denotes the tunnel splitting of ground state. It may be 
worth to estimate the range of validity of our results, i.e., how large m is. QK u sm 2 9o is 
the barrier height of potential and huio is the level space between neighboring levels. For 



the horse-spleen ferritin reported in P2| , P3| the residual spin is S ~ 100 ( corresponding 
moment Mo = 217/ie) and volume is Q ~ 2 x 10 _19 cm 3 (diameter 7.5nm). The longitudinal 
anisotropy constant and transverse susceptibility are seen to be K„ = 2 x 10 5 erg/cm 3 and 
X± = I0~ 5 emu/G cm 3 respectively. Using the above parameters we find that the number 
of the levels in the potential well is about 10 as 5 ~ 0.03. Fig.3(a) shows the oscillation of 
tunnel splittings of lowest 3 states with respect to the external magnetic field due to the 
quantum phase interference of two tunneling paths of opposite windings for Sq =integer and 
half-integer. From Fig.3(a) one can find that the magnitude of tunnel splittings at excited 
states is much higher than that at ground state and may contribute significantly to the 
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experimental observation at finite temperature. When d — So — l — j, i.e., h = (S — l — j)h 
(/ is an integer), the tunneling splitting AE m vanishes. The period of oscillation is 



h 

Ah = — 

I 1 } 



(35) 



which is independent of the energy. 



IV. QUANTUM PHASE INTERFERENCE AS Hy / 



When the external magnetic field is applied along an arbitrary direction in the plane of 
the hard axis and medium axis, the effective potential V(0) — SlK n sin 2 # (sin — b) 2 has the 
asymmetric twin barriers which lead to that Neel vector may rotate from one orientation(0 = 
0+) to another (0 = n— + ) along clockwise underbarrier path and anticlockwise path(Fig.2). 
Two different instantons (Fig.2) corresponding to tunneling through two types of barriers 
are found as 



where 



7T 



±— ± 2 arctan [A±sn(gr, k 2 )] 
2 



{l-e) 2 -b 2 
(l+e) 2 -b 2 



E 



QK tl sm 2 O ' 



(36) 



q =^[{l + e) 2 -b 2 ]\X ± 



[l-e) 2 -b 2 
[l±b) 2 -e 2 



Our starting point for investigation of the tunneling and related quantum phase interference 
at excited states is still the transition amplitude of the barrier penetration projected onto the 
subspace of fixed energy E, i.e., the Eq.(19) from which the tunneling splitting is obtained 

as 



AE 



(2E8\ 



d<t> f d<t>nl}* E {<t>f)i)E{<i>i)G{(j>f, p; (pi, -/3) 



(37) 



The result corresponding to the Eq.(22) is formally the same 
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exD 

ae — *yJ.\i+ + r 2 



(38) 



In the present case , however, the two tunneling paths are not symmetric. Thus we find 



If = exp 



h 



if, 



(39) 



Sf = ±0 [?r =f 2 arcsin(6 ± e)] , 



(40) 



lf = J #±#^(^)^(0f)G' ± (^,/3;0f,-/3), 



(41) 



1 



G ± = I £>[0]exp ( -- J Lf ff dr ) , 



/3 



(42) 



(43) 



2 V dr 

If is now dependent on tunnel direction. The phase dependent tunneling kernel G ± is 
evaluated with the help of the periodic instanton. Following the procedure above we obtain 



/±~2/3exp(-^) 



exp 



Wi- 



ll 



(44) 



a= [(l + £ ) 2 -b 2 ] K k' 2 = v 7 ! - h, 



e x 2 ± 



[X 2 ± E(k 2 ) + (k\ - X 2 ± ))C(k 2 ) - (A 4 ± - k 2 2 )U(k 2 , Xl)} , 



(45) 



where H(k 2 , A±) is the complete elliptic integral of the third kind. Considering the additional 
phase contribution from the real-time paths in potential well the tunnel splitting Eq. (38) 
is written as 



AE = 



If exp 



h 



+ I 2 exp 



iS r 



(46) 



where 
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(47) 



5f = —20 [arcsin(6 ± e) — arcsin b] , 



(48) 



$± = 2/a;o 



E(if± ' k ' 2) + a(l =F 6) 2 F(^, k' 2 ) T 2b*(l T b - £)r%±, a* *£) - e\^\ 
o V l T o 



(49) 



ip = arcsin 4 / -rrz — ,ot = 



2e 



2(1 =F b) \l T b + e' 

Inserting Eqs.(39), (44) and (47) into Eq.(46), we obtain the final formula of the tunnel 
splitting 

Hlo 



AE 



2aJC(k 2 



h 



h 



h 



ra P 1 ~^-) + exp (~^r-) + 2 exp ( + W * ) cos [2A7r - ($+ - $;)] |> (50) 



which is a function of the external magnetic field and the energy. For low lying excited states, 
e << 1 , k' 2 << 1, the energy E is again replaced by harmonic oscillator approximated 
eigenvalues E m — (m + ^)huj. Expanding the complete elliptic integrals E(k 2 ), IC(k 2 ) and 
n(/c 2 , A±) in the Eq.(45) as power series of k' 2 we obtain 



± l<>/\ sin 2 0i) 

UJQ 



^T^¥ - — (1 - b 2 )h 2 4 (\n^r + - ) +6arcsin6^- 
16 \ k 4 / 2 



(51) 



Substituting the Eq.(51) into the Eq.(50) and taking note of limits JC(k' 2 — > 0) — >|, 
a (k' 2 — > 0) — > (1 — 6 2 )^ and ~ $7 at low lying excited states we obtain the tunnel 
splitting of the mth excited state as 

E 2 



AE, r 



ml 



4B(l-b 2 ) 2 [cosh^Evr) + cos(2Att)] 3 , 



(52) 



E 2 = 



ab (i-b 2 y 



exp 



B VT^+b arcsin 



(53) 
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Fig. 3 shows the oscillation of tunnel splitting at low lying excited states with respect to the 
external magnetic field for So =integer and half-integer respectively. When A = (2Z + 1)/2 (I 
is an integer), tunnel splitting AE m tends to a minimum value. The period of oscillation is 

Ah =-*°- (54) 
cos a 

which is independent of the level, but dependent on the direction of the external magnetic 
field. When a = and m = 0, the tunnel splitting AE m reduces to the result in Ref.[15]. 
The period increases with the angle a. When the direction of the magnetic field is along 
the medium axis (a = |), the period approaches to infinity, in other words, the oscillation 
disappears. 

V. CONCLUSION 

The effect of the macroscopic quantum phase interference at excited states is studied 
for the biaxial ferrimagnetic particles with the external magnetic field applied along an 
arbitrary direction in the plane of hard and medium axis. We present a general formula of 
tunnel splitting at excited states as a function of the magnetic field and the energy. The 
oscillation behavior of tunneling splitting at low lying excited states is similar to that in 
FM particles observed experimentally in molecular clusters Fes and should be observed in 
further experiment with ferrimagnetic particles for which a possible candidate of materials 
may be horse-spleen ferritin 
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Figure caption 

Fig.l (a) The periodic potential and the instanton trajectories. The arrow lines denote 
two tunnel paths of opposite windings, (b) The equilibrium orientations of Neel vector in 
the absence of Y-component of the magnetic field. 

Fig. 2 (a) The potential with asymmetric twin-barrier and instanton trajectories, (b) The 
equilibrium orientations of Neel vector in the presence of Y-component of the magnetic field. 

Fig. 3 The level splitting as function of the external magnetic field with angular (a) 
a = 0°, (b) a = 3°,(c)o! = 5° for So =integer(solid line) and S =half-integer(dot line). Here 
S = 100, D, = 10- 19 cm 3 , x ± = 10~ 5 , K u = 10 5 erg/cm 3 and 5 ~ 0.03. 



[1] A. J. Leggett, S. Chakravarty, A. T. Dorsey, M. P. A. Fisher, A. Garg, and W. Zwerger, Rev. 
Mod. Phys. 59, 1 (1987) 

[2] E. M. Chudnovsky and L. Gunther, Phys. Rev. Lett. 60, 661 (1988) 

[3] B. Barbara and E. M. Chudnovsky, Phys. Lett. A 145, 205 (1990) 

[4] L. Gunther and B. Barbara (ed) Quantum Tunneling of Magnetization, QTM 94{K\auwer, 
Dordrecht, Netherland, 1995) 

[5] D. Loss, D. P. DiVincenzo, and G. Grinstein, Phys. Rev. Lett. 69, 3232 (1992) 

[6] E. M. Chudnovsky and D. P. DiVincenzo, Phys. Rev. B 48, 10548 (1993) 

[7] J. -Q. Liang, H. J. W. Miiller-Kirsten, and J.-G. Zhou, Z. Phys. B: Condens Matter 102, 209 
(1993) 

[8] A. Garg, Europhs. Lett. 22, 205 (1993) 
[9] W. Wernsdorfer and R. Sessoli, Science 284, 133 (1999) 
[10] A. Garg, Phys. Rev. B 60, 6750 (1999) 

14 



[11] S. P. Kou, J. Q. Liang, Y. B. Zhang and F. C. Pu, Phys. Rev. B 59, 792 (1999) 

[12] Yan-Hong Jin, Yi-Hang Nie, J.-Q. Liang, Z. D. Chen, W. F. Xie and F. C. Pu, Phys .Rev. B 
62, 3316 (2000) 

[13] J. M. Duan and A. Garg, J. Phys.: Condens. Matter 7, 2171 (1995) 
[14] H. Simanujunlak, J. Phys.: Condens. Matter 6, 2925 (1994) 

[15] Yi-Hang Nie, Yan-Hong Jin, J.-Q. Liang, H. J. W. Miiller-Kirsten, D. K. Park, F. C. Pu, J. 
Phys.: Condens Matter 12, L87 (2000) 

[16] Y.-H. Nie, Y.-B. Zhang, J.-Q. Liang, H. J. W. Miiller-Kirsten, F.-C. Pu , Physica B, 270, 
95(1999) 

[17] E. M. Chudnovsky , J. Magn. Magn. Matter 140-144, 1821 (1995) 
[18] N. S. Manton and T. S. Samols, Phys. Lett. B 207, 179 (1988) 
[19] A. N. Kuznetsov and P. G. Tingakov, Phys. Lett. B 406, 76 (1997) 
[20] J.-Q. Liang, H. J. W. Miiller-Kirsten, Phys. Rev. D 40, 4685 (1992) 
[21] J.-Q. Liang, H. J. W. Miiller-Kirsten, Phys. Rev. D 51, 718 (1995) 

[22] D.D. Awschalom, J.F. Smyth, G.Grinstein, D.P. DiVincenzo, and D. Loss, Phys. Rev. Lett. 
68, 3092 (1992) 

[23] J. Tejada and X. X. Zhang, J. Phys.: Condens. Metter. 6, 263 (1994) 



15 




fig.2 




(b) 



Fig.3 




Fig.3 



